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It has been recently proposed that a sterile neutrino of mass = 40-80 MeV, 
mixing ~ 10~ 3 -10~ 2 , lifetime Ss 10~ 9 s, and a dominant decay mode — > 

could be the origin of the experimental anomalies observed at LSND, KARMEN and 
^ \ MiniBooNE. Such a particle would be abundant inside air showers, as it can be pro- 

duced in kaon decays (K — > Vhfi, Kl — > v^t:^). We use the Z- moment method to 
evaluate its atmospheric flux and the frequency of its decays inside neutrino telescopes. 
We show that would imply around 10 4 contained showers per year inside a 0.03 km 3 
telescope like ANTARES or the DeepCore in IceCube. These events would have a 
characteristic energy and zenith-angle distribution (E u ~ 0.1-10 TeV and 9 < 90°), 
which results from a balance between the reach of the heavy neutrino (that disfavors 
low energies) and a sizeable production rate and decay probability. The standard back- 
ground from contained neutrino events (v e N — > eX and neutral-current interactions of 
high inelasticity) is 100 times smaller. Therefore, although it may be challenging from 
an experimental point of view, a search at ANTARES and IceCube could confirm this 
heavy-neutrino possibility. 



1 Introduction 



The direct observation of neutrino interactions in different types of experiments pQ has been 
used to establish a basic picture of neutrino masses and mixings. From a model building 
point of view, this is arguably the most significant discovery occurred in particle physics 
since the confirmation of the standard model in the early 70 's, as it reveals a scale that is 
(most likely) not electroweak. The picture, however, has faced some persistent anomalies 
in experiments with neutrino beams from particle accelerators. Basically, muon neutrinos 
of energy below 1 GeV seem to experience an excess of charged-current (CC) interactions 
with an electron in the final state. The interpretation of these events in terms of — > v e 
oscillations is inconsistent with the mass parameters deduced from solar, atmospheric and 
reactor neutrino observations. 

In a recent analysis Gninenko [2] has made a very compelling case for a massive neutrino 
as the origin for all these anomalies: 

(i) LSND [3] observed i> e -like events (z? e P — >* e + n) with a gamma signal from neutron 
capture that seem to imply an excess of — > v e oscillations. He shows that the events could 
be equally explained through production (v} 2 C — » vyjiX) followed by its radiative decay 
Vh — > v li with the final 7 converted into a e + e" pair indistinguishable from an electron. This 
explanation would work for a large enough production cross section (\U^h\ 2 ~ 10~ 3 -10~ 2 and 
m u < 80 MeV) and a short enough decay length (r^ ^ 10~ 9 s and m v > 40 MeV). During 
the first years of data taking LSND also observed an excess of v} 2 C — > e~X events that 
were interpreted as v ^ — > v e oscillations but are consistent as well with the hypothesis. 

(ii) KARMEN [1], using a similar technique, did not confirm the LSND anomalies. The 
neutrinos at LSND, however, had an average energy of 100 MeV and a long high-energy tail, 
whereas the spectrum at KARMEN was a narrow peak around 20 MeV. Gninenko shows 
that a 40 MeV neutrino would be above the production threshold at KARMEN, which makes 
his hypothesis consistent with the data. 

(in) MiniBooNE [5] has observed an excess of electron-like events for v e energies between 
200 and 475 MeV, with no significant excess at higher energies. Gninenko's fit exhibits also 
a good agreement with the data (higher-energy events are disfavored by an increase in the 
decay length and are hidden by the low statistics). More recently [B] this experiment has also 
reported an excess in data for antineutrino energies in a wider range. His fit is consistent 
as well, and could favor a Dirac nature for Vh. 

In addition, the mass range 40 < rrih < 80 MeV makes Uh too heavy to be produced in pion 
decays and too light to distort the muon spectrum in kaon decays. The heavy neutrino is also 
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produced when the muon itself decays, but he argues that a mixing \U^h\ 2 < 10~ 2 makes it 
acceptable. Specific searches for unstable neutrinos put strong constraints on {U^l, but are 
based on decays with charged particles in the final state (u h — > eeu, /xeu, (jlttu), never on the 
decay Vh —> ^7 induced by a magnetic moment transition. Its large mass and short lifetime 
should keep Uh also safe from bounds from supernovae and primordial nucleosynthesis [7]. 
Finally, a recent analysis [5] of muon capture with photon emission at TRIUMF finds that 
Gninenko's neutrino would imply a signal well above the 30% excess (versus the standard 
model value) deduced from the data One should notice, however, that the photon energy 
cut and the small size of the target volume at TRIUMF make this experiment very sensitive 
to the neutrino lifetime. A value Th ~ 3 x 1CT 9 s could imply a consistent radiative capture 
rate there while explaining the data at LSND and KARMEN (which require < s) 
and still having an impact at MiniBooNE0. 

We find the heavy neutrino hypothesis very interesting and will study here its implications 
in a different type of experiments. Our basic observation is that would be abundantly 
produced in the atmosphere through kaon decays. At energies around 1 TeV its decay 
length becomes cr^7 rj 5 km, which implies that can reach a neutrino telescope and 
then decay. The final photon would be seen there as a 'pointlike event, similar to the shower 
from v e N — > eX or from a neutral-current (NC) interaction of high inelasticity but clearly 
distinguishable from the muon track in v^N — > fiX. 



2 Neutrino fluxes at sea level 

The atmospheric flux of any species can be easily estimated using the Z-moment method 
[TOT ITT] . This method provides a set of coupled differential equations that describe the 
evolution with the atmospheric depth t (in g/cm 2 ) of the fluxes of parent hadrons ((f) h with 
H — p, n, tt^, K^, Kl) and of any particles that may result from their decay or their collision 
with an air nucleus. The generic equations for (f>H(E,8,t) are 



- Xjf + $H'H , (1) 



^ -^dec ^int H 

where A^ c (A^ t ) is the decay (interaction) length of H in the air and Sh'h describe the 
sources. These equations can be solved analytically under some simplifying assumptions, 
namely, 

(i) the all-nucleon primary flux has a constant spectral index —a; 



X A global fit including TRIUMF would certainly constrain further the parameter space in Gninenko's 
model. 
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(ii) the energy distribution of particles from collisions and decays scales linearly with the 
energy of the parent hadron; 

(in) the hadronic interaction lengths A^ t do not change with the energy; 

(iv) the contributions to the nucleon flux from meson collisions and to the pion flux from 

kaon collisions are negligible. 

It follows that the nucleon fluxes 0tv keep the same spectral index —a at any depth, and 
that the source terms are reduced to 

Snh = tjt Znh (2) 

A int 



dx x a - x F NH (3) 



where the Z-factors 

7 / J™ ™a— 1 

^>NH — I 
JO 

are constants (independent of E and the zenith angle 9) derived from the distribution F^h(x) 
of the fraction of energy taken by H after a iV-air collision. The meson fluxes 4>m can then be 
easily derived in two different regimes. At low energies A^ c = (E/rriM)cTMP is much smaller 
than A^f t and meson interactions can be ignored, whereas at high energies the variations in 
<Pm are dominated by collisions with air nuclei (the air density p is a function of t and of 
6). A simple interpolation can be used between these regimes. We take in our analysis the 
primary flux, the Z-factors and the atmospheric model in [TT] (see [32] for a discussion of 
the fluxes at higher energies). We obtain, for example, that the TeV charged-pion vertical 
flux reaches its maximum (a 4% of the initial nucleon flux) at a depth of 200 g/cm 2 , and 
that the kaon flux there is 7 times smaller. 

The lepton fluxes from meson decays can also be incorporated. In particular, standard 
neutrinos do not interact nor decay in the atmosphere and their flux equations will only 
depend on source terms of type 

S Mv (E) = B(M v) \\xx- 1 *^J*\ f Mv (x) , (4) 

where B(M — y v) is the branching ratio of a given decay mode, Fmv{%) is again the distri- 
bution of the fraction of energy taken by the neutrino in that decay, and the dependence 
on t and 9 is implicit. We obtain that, although kaons are less abundant than pions in air 
showers, a lower ratio A^ c / A^ t makes them the main source of neutrinos at energies above 
100 GeV. The TeV flux at sea level is dominated by muon neutrinos, with <j) p w 0.42 0^ , 
<j) Ue w 0.036 0^ and Pe w 0.023 0^. 

Heavy neutrinos will be mainly produced in charged-kaon decays. The branching ratio 



is 



B (K + -> p + v h ) w B{K + -> p+v) x |f/^| 2 p h , (5) 
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Figure 1: Neutrino fluxes {i>i + z/j) at sea level for 9 = (solid) and 9 = 60° (dashes). 

where B(K + — > p^v) = 0.64 and the kinematic factor for m^ = 40-80 MeV is ph ~ 
(1 + fn\/m 2 ^) [13]. The fraction of energy x taken by will have a flat distribution (a 
constant F Kh (x)) between x min and x max , 
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There will be smaller contributions from K + — > 7i°fi + Uh and Kl — > ^ p + Vhi plus the analo- 
gous K~ and Ki decays into Vh (the heavy neutrino may be a Dirac or a Majorana particle, 
see discussion in [2]). The equation defining (f> h (E,9,t) is 



dt 



A dcc 



Kh 



(7) 



A 



where the source terms take the form in Eq. (0J, A^ ec = (E/mh)cThp, and the sum runs over 
the decay modes that produce v h . In Fig. 1 we plot the total heavy neutrino flux at sea 
level from inclinations 9 = 0°, 60°. We have taken the central values m^ = 60 MeV and 
|[/ M | 2 = 0.005, with r h = 10" 9 s. At 1 TeV 73% of the flux comes from K + decays, K~ 
contributes a 25%, and K L just a 2%. Finally, notice that the photons produced in the air 
through Vh decays are together with other photons and muons inside the parent shower and 
are therefore non-observable. 



5 



Figure 2: The total flux through the lateral surface of the detection region cancels if Rt ^> 
dj cos 6 (we take 9 < 85°). 

3 Events at a neutrino telescope 

As neutrinos enter the ground their sources disappear and they just experience two types 
of processes: heavy neutrinos may decay into 71^, whereas and v e may have neutral- 
or charged-current interactions with matter. At a depth d the sea-level fluxes <f)i(E,8,0) 
become 



where d and the decay/interaction lengths are given in meters and we have neglected the 
curvature of the Earth (a good approximation for 6 < 85°). At 1 TeV we obtain \\ ec ~ 5 
km whereas Aq C ^2x 10 6 km and A^ c »8x 10 6 km (the interaction lengths decrease with 
the energy as I/c^tv). This means that the decay of a Uh crossing a neutrino telescope is 
10 6 times more probable than the interaction of a standard neutrino. In addition, CC 
interactions will be clearly different from decays, as the final muon will produce a track 
hundreds of meters long [15] . The electromagnetic shower from a Vh decay will be pointlike 
(it develops in a few meters), similar to the one produced by a NC interaction or a v e CC 
process. These standard events, however, are suppressed by the lower u e fluxes [y^ — > u e 
oscillations at L < 100 km and E > 100 GeV are negligible) and the inelasticity distribution 
(oc 1/y [14J) in v-N collisions. 

Let us be more specific. To estimate the number of events per unit time occurred inside 
a telescope one needs to calculate the total flux (ingoing plus outgoing neutrinos) through 
the surface containing the detectors. We model this region as a cylinder of section A and 
length H starting at a depth d (i.e., d goes from d Q to d + H). For a fixed angle 6 ^ 85° 
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Figure 3: Energy distribution of the contained showers from — > 1/^7, uN — > vX and 
v e N — > eX at ANTARES. We also plot (dashes) the energy distribution of the parent 
neutrino in each case. 

the neutrino flux only depends on the depth d. Therefore, the total flux through the lateral 
surface of the detector will be zero (given 9, the flux through any lateral dS is equal to the 
flux leaving the detector through an opposed lateral surface — dS, see Fig. 2). The number 
of heavy-neutrino events inside the detector in an interval of energy and solid angle per unit 
time can then be calculated as the difference between the fluxes through its upper and its 
lower sections: 

N h = [ dE ( dtt ( dS cos 0U h (E,G, do) 

J AE J AQ J A 

ME,0,d o + H)]. (9) 

An analogous expression can be obtained to estimate the number N u of interactions inside 
the detector produced by neutrinos from those directions. In Fig. 3 we plot in dashes the 
energy distribution of the neutrinos that interact (v^ and v e ) or decay (z/^) per year inside 
a detector like ANTARES p] (A = 0.1 km 2 , d = 2.2 km, H = 0.3 km) or the DeepCore 
[T7] in IceCube (of similar size and depth). The energy of the initial neutrino, however, is 
not the most relevant parameter for observation, as in NC interactions only a small fraction 
y may be deposited in the detector. In a vh decay only the photon energy is visible (we will 
assume an isotropic decay [2]), whereas in i/ e -CC interactions all the energy carried by the 
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neutrino goes to the contained shower. If the inelasticity y in the event has a distribution 
F us h(y), then the energy distribution of the cascades inside the detector will be 



Fig. 3 shows in solid lines the main result of our analysis. The number of standard showers 
of energy above 100 GeV from down-going neutrinos inside ANTARES is around 1300 per 
year (60% from z/ e -CC interactions and 40% from NC interactions). For the central values 
m h = 60 MeV and {U^ 2 = 0.005 the heavy neutrino would provide 26000 extra events. If 
the energy threshold is set at 500 GeV the number of standard events is reduced to 220 per 
year, whereas the number of events from decays is just cut to 14000. 

4 Summary and discussion 

Telescopes like ANTARES or IceCube are designed to observe upward-moving muons pro- 
duced in neutrino interactions near the detector. These events are clean, in the sense that 
no particles except for neutrinos can reach the detector after crossing the Earth. Telescopes 
can also observe the contained showers produced in NC interactions or in z/ e -CC processes. 
Since their development takes just a few meters, these events are pointlike, and the only sign 
indicating whether they are caused by an upward or a downward-going neutrino is that in 
the latter case they may come together with muons. 

In this paper we have shown that the decay of a long-lived neutral particle produced 
in the atmosphere could change drastically (by over a factor of 100) the number of TeV 
contained showers in these experiments. In particular, we have analyzed Gninenko's heavy 
neutrino, that appears as a possibility well motivated by the results at LSND, KARMEN and 
MiniBooNE. We find remarkable that its mass, mixing, and lifetime optimize the distortion 
introduced in TeV-neutrino telescopes: below ps 100 GeV Vh does not reach the telescope, 
and above 100 TeV its decay length becomes too large and the signal vanishes (notice that 
A^ ec grows with the energy while A^ nt decreases). 

The heavy neutrino would be produced through kaon decays inside air showers together 
with a muon of similar energy. A crucial question is then whether these decays can be 
disentangled from the muon bundle associated to the parent shower. ANTARES or the 
DeepCore in IceCube are more than 2 km deep, and as the zenith angle grows all muon 
effects will decrease. In contrast, the zenith angle dependence of the vh events up to 85° is 
very mild (especially at larger energies), since these neutrinos do not lose energy in their way 
to the detector. Therefore, the TeV contained showers should appear as a clear anomaly that 




(10) 
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may be accompanied by muons in vertical events but that persists at higher zenith angles 
(slant depths). The large number of events that we obtain could allow for specific searches. 
For example, events with lower-energy muons in the upper part of the detector followed by 
a contained TeV shower below, or other topologies that would otherwise be discarded. 

If a heavy neutrino is the explanation of the LSND and MiniBooNE anomalies, our results 
show that it will reach effectively the core of neutrino telescopes and will decay there at a high 
rate. A MonteCarlo simulation of individual showers, including all the muon backgrounds 
and the response of the detector, should then provide the best strategy in the search for an 
observable signal that confirms or excludes this heavy neutrino possibility. 
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